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Motivation 


Fully automated meshing for Reynolds-Averaged Navier-Stokes Simulations 


Mesh generation for complex geometry continues to be the 
biggest bottleneck in the RANS simulation process 


Embedded boundary Cartesian methods routinely used for 
inviscid simulations about arbitrarily complex geometry 


These methods lack of an obvious & robust way to achieve 
near wall anisotropy 


Goal: Extend these methods for RANS simulation without 
sacrificing automation, at an affordable cost 
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* Previous work & analytic wall functions 


* ODE-based wall models 


- A New ODE wall model 


* Numerical examples 


* Conclusions 


Previous Work 


Analytic wall functions 


* Thin-layer form of streamwise a) ( =) Op Ou n | 


momentum for RANS eqs. (ut M5, Or +p it 
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Previous Work 


Analytic wall functions 


+ Thin-layer form of streamwise G) Ou Op Ou Ou 
(e+ Het) 


momentum for RANS egs. Oy ay an +? |\“de t "Oy 


+ The diffusion model assumes that 
velocity is small and ZPG 
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Previous Work 


Analytic wall functions 


Thin-layer form of streamwise G) Ou Op Ou Ou 
iS : (M+ He) 3 = +pluz—+v; 
momentum for RANS eqs. Oy ay ar On oy 


The diffusion model assumes that Velocity profile 


velocity is small and ZPG 30 
— Turbulent flow 
9 — SA wall function Outer. 
é Ou 25{ _— Spalding profile flow 
5 | (ut ui)5— ] = 0 
oy Oy 
20 
+ Assume mixing-length model for eddy is 
£15 
viscosity: vt ~ distance to the wall ia 
Le = pur = pkvyt 10 
* Gives a very good fit to experimental 5 
velocity data up through the log layer 
10° 10" 10" 0° 10° 


Previous Work 


Analytic wall functions 


Spalding model: 
> : 1 : 1 
yt (ut) =ut + e788 Gal —1— «ut 5 (nu y= ; (Kut 5) 
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SA wall function (2012): 


Derived, using a limiting form of SA turbulence model and integrating the diffusion model 


Previous Work 


Analytic wall functions 
Velocity profile 


* Spalding model: yt (ut )=... — Turbulent flow 
— SA wall function 
= 25| .—— Spalding profile 


+ SA wall function (2012): ut (yt) 


+ Both: a) 


Wake 


+ Are good approximations and give accurate 2 region 
wall shear stress when anchored with F’ 215 Log | 
ut(yt) = B+ ec log((yt +41)? +7) — c2log((y* + a2)? + b3) located out to the log-layer faver 
10 
_ so Hata Ohare =. fy AvetanD ye dee Be + Are inappropriate beyond the log layer (in 
c3 arctan2(y* + ay, b1) caarctan2(yT + ae, b2) the wake region) : 
5 Transition 
+ Prefer SA wall function, since it gives direct relationship for velocity as a function of distance region 
* Knowing wat a point F, iterate to find ur, so that ut(yf) = up = urur 10° 10! 10" 10° 10° 
Previous Work Previous Work 
Constructing forcing points Analytic wall functions 
30 008 = Fine er 
. i ima — SA wall model Sm ee 
\ | + 7 2s{ _— Spalding profile 0.0 
=e 2 2 p |. yt | _ ooo 
F 20 


* Construct forcing points at uniform distance from wall 
+ Interpolate data to point F from cell centered solution on outer grid 


* With velocity an distance at forcing point, use wall function to find ur 
and wall shear Tywail = pu 
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Pressure Coefficient 
Skin Friction Coefficient 
5 


— Fine 
© Experiment (Cook et) 


rs 
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‘yor 0s “03 0a cr 06 07 08 09 1 
* Successfully coupled applied analytic wall functions with cut-cell Cartesian meshes in 2012* 
* Introduced new SA wall function - worked well where wall functions were appropriate 
* Results were comparable to body-fitted methods using wall functions 
* Conclusion: 
Cartesian RANS is viable, but wall functions alone are probably not 
sufficient to make the approach cost competitive 


* AIAA 2012-1301, "Progress Towards a Cartesian Cut-Cell Method for Viscous Compressible Flow", Berger, Aftosmis & Allmaras 


Good wall functions gone bad 


==" Mixing-length model 
—— Eddy viscosity 


— «velocity 


Thick boundary layer 
. UP = 10? 
* Forcing point in log layer 


* Mixing length model gives good 
estimate eddy viscosity 


¢ Analytic wall function is appropriate 


Good wall functions gone bad 


==" Mining-lengh modst 
Eddy viscosity 
velocity 
Thick boundary layer 

. UP =~ 10° 

* Forcing point in log layer 


* Mixing length model gives good 
estimate eddy viscosity 


¢ Analytic wall function is appropriate 


Thin boundary layer 


+ At the same distance from the wall, UPR ~ 104 
* Forcing point is now in the wake layer 
* Eddy viscosity highly non-linear 


* Mixing length model is a poor approximation, analytic wall function 
inappropriate 


Outline 


* ODE-based wall models 
- A New ODE wall model 


ODE-Based Wall Models 


Proposed by several authors in last decade 


. Tangential 
* Solve ODE on 1D “linelet” normal to surface UaUSA Teo 
* Solve: 
. . ‘ Eddy 
* Diffusion eq. for streamwise momentum viscosity 


* Turbulence model in wall-normal direction 


* Produces a system of 2-point, 2nd-order 
BVPs 


1D linelet 


N y 
Wall 
* Coupling: Just like an analytic wall function 


see Kalitzin et al., J. Comp. Phys., 204, 2005, Bond & Blottner, Intl. J. Num. Methods Fluids, 66, 2011, or Capizzano, AIAA J. 54(2), 2016 


ODE-Based Wall Models 


SA-BVP: Diffusion equation coupled with wall-normal SA turbulence model 


«-momentum: fa) 
(diffusion equation) 


SA model on linelet: 


a _. Ov av\* 
- (v +0) oe = -Cp2 oe + Production — Destruction 
Oy Oy Oy 
ee 
— 
wall-normal 


diffusion 


ODE-Based Wall Models 


Compare SA-BVP with SA wall function on turbulent bump in channel 


Mo = 0.2 
Ret = 3x10° 


a) Ou 
a@-momentum:  =— { (w+ pe) 
Oy Yy 


a » 8lasn *) 
A model on linelet: By F ay 


ay 
Ov\~ 
—Cp2 ( ) + Production — Destruction 


ODE-Based Wall Models 


Compare SA-BVP with SA wall function on turbulent bump in channel @ z = 1.2 


Forcing 
int 


whtedge 


— CFL3D (fine grid) 
=~ Mixing length model 
— SA-BVP: eq.(11) 


— CFL3D (fine grid) 
~~~ Mixing length model 
— SA-BVP: eq.(11) 


0 omes oot ams os % ams oT 
Distance from wall 


oms “ona ines 
Distance from wall 


* Forcing point well out in wake layer, y= 0.012, w= 0.85 Uedge 


« Mixing length eddy viscosity inappropriate, so diffusion model alone does poorly 
* Improved eddy viscosity makes a significant difference 


ODE-Based Wall Models 


Streamwise momentum in the wake 


Boundary layer profiles @ a= 1.2 


* 
0.01 
Mo = 0.2 
Rey = 3x108 a 
0.0) 
—* 
— = 
~ = 
freaeeeente + 0.01 7m 
I t t | 7 Pe 
ty) 0.4 0.8 1.2 16 a 
0.02 Pr/p+ uur + vuy 1 
te) 0.01 0.02 0.03 


Distance from wall 


* Thin layer streamwise momentum: ((/1 4 Hi )Uy),y Dx + p(t, + Vy) 


« Examine relative magnitude of terms as we move away from the wall 


ODE-Based Wall Models 


Streamwise momentum in the wake pies 


Boundary layer profiles @ a= 1.2 


« At the wall, we have 
(te + fr )Uy), = Pa 


0.01 


* Outside the boundary layer we approach: 90 


0 = py + p(uttg + vty) 


+ In between we have the full streamwise 
momentum eq. 


u/10 
ur 
vay 
Pip + uur + vay 


Pa + p(urty, + vty) ~o 
LI LI LU 


(Cu + Hey), 


0.01 0.02 
Distance from wall 


0.03 


ODE-Based Wall Models 


Streamwise momentum in the wake oes 


Boundary layer profiles @ a= 1.2 


Forcing point at y = 0.012 


ak Forcing oot 


point : 
3,06 i 
3 0.0 : 
3 : 
FS 
oa ’ 
— CFL3D (fine grid) E 
=~ Mixing length model on * ——= Wid 
02 : — ?, 
— SA-BVP: eq.(11) i it 
; — ty 
aos 0 + ree wee 
- i 10) 0.01 0.02 0.03 


Distance from wall 


Distance from wall 
* Forcing point @ y = 0.012 is in the wake. (4+ 1)Uy), = Px + 6 
¥ x 


The convective balance has similar wu 


magnitude as pz — Need to include! 
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ODE-Based Wall Models 


Including the convective balance 


* At the wall, velocity is zero, convective balance is zero 


(0+ He)uy), = Pe 


* But at the forcing point, it has the same magnitude as pz 


(H+ pe)Uy)y = Po + p(uter + Vty)| jp 


* Computing wall-normal variation of convective balance 
introduces streamwise coupling, and means computing 


the wall-normal velocity — prefer not to do this 


Tangential 


U, VSA 
locity 


Eddy 
viscosity’ 


1Diinelet 


Wall 
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ODE-Based Wall Models 


Including the convective balance 


* Introduce a cutoff function, wy), to turn it off as we 
approach the wall 


0 (as you _ op , »(y) oul _ oul 
Oy aac Oy) Aa" PUIP VUP a rl a | °F Oy FP 


a] ~, OV ices : : 
=—|(vt+v) = = wall-normal diffusion + Production — Destruction 
oy Oy 


where tp(0) = 0, (yr) = 1 


A New ODE-based Wall Model 


bup4 wall model 


* Introduce a cutoff function, y( y)s to turn it off as 
we approach the wall 
usa(y) = 


ugaly) 


USAF 


USAF 


sothat (0) = 0, and (yr) = 1 


Tangential 
locity 


=. sa 


U, VSA 


Eddy 
viscosity’ 


1Diinelet 


Wall 
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A New ODE-based Wall Model 


bup4 wall model 


U, VSA —_ 
Ki vel locity - 
* Introduce a cutoff function, wy). to turn it off as ._ 
we approach the wall 
usa(y)  Ugaly) viscont) 
w(y) = USAW sA\l 
USAF USAR 


1Dlinelet pp 


sothat (0) = 0, and (yr) = 1 


Wall 


* The complete bup4 model becomes 


0 ( you _ Op + Ou cae Ou 
Oy eke Oy) Ox * ve ural ale 


~, OW 
z ( +7) mz) = wall-normal diffusion + Production — Destruction 
y y 
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A New ODE-based Wall Model 


bup4 wall model: Include streamwise convective balance and pressure gradient 


point id Forcing 


point 


7 
7 
250] t Forcing 
i 
i 
i 


2 
3 406 
2 130 3 
ia = 
Ss = 
3 oa 
100) 
— CFL3D (fine grid) —— CFL3D (fine grid) 
=~ Mixing length model ial ~~~ Mixing length model 
baad — SA-BVP: eq(11) — SA-BVP: eq(11) 
% 0.005, ‘001 ‘0015 O02 0025 ay 0.005 oor ‘01S ‘002 0025 


Distance from wall Distance from wall 


on a&—momentum: p (w+ a ap +y(y)p |r oul + | 
= > —((utm)—)=— +0 Up — fala 
ay (eT ay) = ae TPP Parle TOF Gale 
; . Od ~, OV igen P : 
SA model on linelet: iy (v+D) By = wall-normal diffusion + Production — Destruction 
y y 
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Eddy Viscosity, vr 


A New ODE-based Wall Model 


bup4 wall model: Include streamwise convective balance and pressure gradient 


Forcing 
point 


— CFL3D (fine grid) 
~~~ Mixing length model 
— SA-BVP: eq.(11) 
— BVP4: eq.(14) 


— CFLS3D (fine grid) 

= =~ Mixing length model 
— SA-BVP: eq.(11) 
— BVP4: eq.(14) 


boos om amis 902 
Distance from wall 


0 oms oot rs os 
Distance from wall 


bup4 Op 


az-momentum: 2 (ut jou + adly)o|e =| “4 ou 
~ +o MET Ht) = p Up eas 
oy ee oy Ox Whee Youle” OylF 


SA model on linelet: — ( +0) x) = wall-normal diffusion + Production — Destruction 
y 
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A New ODE-based Wall Model 


ODE solver for bup4 


(2) ie sou Op v(y) eal ou 
= —)=—+y pales ees 
ay (+ pe ay Bn + (y)p | ur Bale + Up aylr 


) ~\ OV Pee ' : 
= ((v+7) — } = wall-normal diffusion + Production — Destruction 
Oy Oy 


* Reformulate 2nd order equations as system of four 1st order BVPs 
* Solve with 6th order adaptive ODE solver from Shampine and Muir 
« Use warm starts on each linelet after initial solve ~2x cost of analytic WF 


* Other details of implementation and coupling in paper 
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Numerical Results 


Verification and Validation using examples from the NASA Turbulence Modeling Resource 


* Computational Examples from TMR 


1. Turbulent bump in channel 
2. NACA 0012 


3. NACA 4412 with trailing edge separation 
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Turbulent Bump In Channel 


TMR: “VERIF/2DB: 2D Bump-in-channel Verification Case” 


Mo = 0.2 / 
Ret = 3x 10° | 
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Turbulent Bump In Channel 


TMR: “VERIF/2DB: 2D Bump-in-channel Verification Case” 


Mo = 0.2 { 
Rey = 3x 10° | 
z= 0.2 r= 0.75 | f= 18 
t t t t {> £ 
0 0.4 08 12 1.6 


* Inlet & exit 25 units away, symmetry plane 5 units above 
* Mesh-converged body-fitted results on 1409 x 641 mesh (~900k points) 


* Compare results with CFL3D reference solution with SA turbulence model on 
finest mesh 


Turbulent Bump In Channel 


Bump: lsobars and surface pressure comparison 


Cc 
P98 05 04-05-02 0.1 0 01 02 03 sige — CFL3D 
—— 11 level mesh 


06r 


“02 04 06 O08 T Zz 14 16 
X-Distance 


# Cells or Points Minimum Cell Dimensions 
Coarse 10 levels 52k Ay = 0.0012 
Medium 11 levels 70k Ay = 0.0006 
Fine 12 levels 91k Ay = 0.0003 
Reference solution 900 k Ay =2.5e-7 


Turbulent Bump In Channel 


Bump: Eddy viscosity, vsa VSA 
0 20 40 60 80 100 120 140 160 180 200 


ODE Wall 
Model 


CFL3D 
reference 
solution 


0 05 1 18 
x 


* Good agreement for evolution and peak eddy viscosity 


* Slight negative values of vsq outside of boundary-layer due to 2nd-order 
advective terms, easily controlled by negative-SA turbulence model 


Turbulent Bump In Channel 


Bump: - Skin friction & y* distribution comparison with bvp4 


0.008 ——CELSD (1409 x 641) 
—— lev (coarse) — 1kevel 
— tev (medium) 
Nor, 12 lev (fine) 400 


—— Illevel 
—— 12kevel 


0 025 «205 ~~ «(0.75 I 125° 15 0 025 05 — 0.75 I 125 (15 
X-distance X-distance 


* Smooth Cf historically challenging for cut-cell meshes, but look good here 
* Slight noise from HLLC flux when face-normal velocity passes through zero 


* Good agreement progressing toward mesh convergence, results ordered by 
dissipation 


Turbulent Bump In Channel 


Bump: Boundary layer velocity and eddy viscosity profiles 


g= U2 f= 12 
Mo = 0.2 —— /_, 
Ret = 3x10® = —— ——, 
SSSSeesse SSS 
| [> 2 
0 0.4 0.8 1.2 1.6 


* Compare profiles at three stations 
* Boundary layer thickens by approx. factor of 2 at each station 


* Since resolution of Cartesian mesh is constant, resolution roughly 
doubles each time we move downstream 


ee ee ee ee 


Turbulent Bump In Channel 


Bump: Boundary layer velocity and eddy viscosity profiles 


Tangential velocity Eddy Viscosity 
y outer y outer 
= t 2 Cartesian Cartesian 
mesh ——, 
| 
_—________ 5] 
H__ >} 
i / 
i—_—_»>_ 


tt 


1.2 1.6 


Turbulent Bump In Channel 


Bump: Boundary layer velocity and eddy viscosity profiles 


Next comes SA vs BVP4 at 
coarsest station 


Then comes coarse/medium/ 
fine -- but need to mark pt F 


Turbulent Bump In Channel 


Bump: Boundary layer velocity and eddy viscosity profiles 


Turbulent Bump In Channel 


Bump: Boundary layer velocity and eddy viscosity profiles 


02 02 
0 025 025 
CFLAD (1409 x 641) — cFL3D (1409 641) 
2 Tote (sane) 2 2 TS Gaetan z 
a ialevtim) ane pw ce = taki 3° 
z Walrmoictaoey |S z 3 Wallmodsi (ios) | 3 
3 ~ Wall model (11 lev) 3 2 doas ‘Wall made (1 ev) Zous 
= ---- Wall model (12 lev) 2 > x 2 Wall model (12 lev) cf 
E01 £01 3” i 5 
€ & a Eo Eo 
wo x= 0.75 
00s oos = 
oos| zo 005 
0 2 oF ry = + > z a i> e 2 oF @ = 7 7 = z 2 
: OO vcrance hom val ue ” m 10 ancoeciawel: on av ota Sa vstanee toms val ee " MO staiies Pecnieil wv : Ot Guechoned 
* 7 curves on each plot, (wall model & field solution) x (coarse, med, fine) + CFL3D * About twice as much Cartesian resolution at a = 0.75 
* Very good agreement for velocity, good agreement for eddy viscosity * Profile shows effects of moderate favorable pressure gradient on the front of the bump 
* z= 0.2 is the most under resolved station, ~4-5 Cartesian cells in boundary layer * Data from the wall model collapses to reference solution regardless of outer resolution 
¢ Eddy viscosity slightly overpredicted due to lack of resolution in outer mesh 
a 28 
Turbulent Bump In Channel S Turbulent Bump In Channel 
Bump: Boundary layer velocity and eddy viscosity profiles Bump: Compare bvp4 with analytic SA wall-function 
Skin friction on coarse mesh 
02 02 ube cats a09 out) 
i BX Van fanetion 
— CFL3D (1409 x 641) baer 
— 10 lev (coarse) 
Boas = ites (medium 20.5 0006 
q 12 ev (ine ; 
3 ite gs |: iia 
S Wail model iiiew) JS e 7 
Zo, == Wall model (12 lev) Zo, Sonne 
&, &, 
FS 2 0.003 
os e=1.2 os 0.002 
5 b> 
renee eeeeterer ono Re 
F F : Bates [eS 
a rT a aT aT rT a ee aC a a eC a or eas os Toys ras Ts 
Distance from wall Distance from wall Distance from wall ctistanee: 
* Aft of bump, slight adverse pressure gradient, thick boundary layer + Even on coarsest grid, SA wall function does reasonably good job 
* " A i ~~ apt mw 
* Velocity profiles show very good agreement -- even on semi log scale * To see differences, look up front on coarse grid, = 0.2, (y* = 280) 
* Eddy viscosity peak being eroded slightly by dissipation on outer mesh 
20 40 
4 
H 
Turbulent Bump In Channel Turbulent Bump In Channel H 
: mixing-length 
Bump: Compare bvp4 with analytic SA wall-function Bump: Compare bvp4 with analytic SA wall-function {model for eddy- 
: viscosity 
Skin friction on coarse mesh Skin friction on coarse mesh ’ 
‘tie cep hws) Salama ca ab 09 xi) 0.008 i oats xem) — crap (1400 x 641) : — Sun any 
= SA Wal funtion 02[- > $A wal funtion os Sh pall tion eee oa] == SA wal ct : Ss SA al neton 
= — > Wall Model a s=-0 "Wall Model: = — > Wall Model “ ‘eo Wall Model. 
zB B 
Boas 230 301s Bu 
s E s E 
~ z 2 > z Z 
é a é a 
0 0 
00s 00s 
: x= 0.2 of ¢ = 0.2 : x= 0.2 of c= 0.2 
“TS for oF to ao? to? 0? oO “Ont O1902‘TNs” O-00F 0905. 01906 0.007 008 “TS {7 a0? ao ao? a0? ODOT -0.003-'0.003- 0.003 0:905~ 0.006 -0.007-_O.08 
X-distance Diitence fom wall Distance from wall X-distance Distance from wall Distance from wall 


« Even on coarsest grid, SA wall function does reasonably good job 
* To see differences, look up front on coarse grid, «= 0.2, (y* = 280) 


* Skin friction discrepancy comes from misprediction of eddy viscosity by analytic wall 
function since it assumes a mixing-length model 
at 


* Even on coarsest grid, SA wall function does pretty good job 
* To see differences, look up front on coarse grid, 2 = 0.2, (y* = 280) 


* Analytic wall function overpredicts eddy viscosity by about factor of 3, 
—is inconsistent with outer solution 
42 


NACA 0012 


Modified NACA 0012 geometry with sharp trailing edge 


NACA 0012 


Modified NACA 0012 - Surface pressure mesh convergence 


— CFL3D (14.7M cells) 


—= 16 level (medium) 
— 17 level (fine) 
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* 1. NOW COMES vs SA-BVP 
* (slides 49/50, but first show skin friction, then show magnitude of terms) 
« show what different levels of modeling are bying us 


* 2. Then comes the cut at 1.001 
* 3. Move adaptive to backup 


6 — 15 level (coarse) 
— 16 level (medium) 
Me = 0.15 Mo = 0.15 — 17 level (fine) 
Rey = 6x10° Rew = 6x10° 
Geo — 10° Qo = 10° 
a = 
— he 
— Bes 
— — 
— 
0 32 oF a6 os 
x/C 
* Validation example “2DN00: 2D NACA 0012 Airfoil Validation Case” of TMR website # Cells or Points Minimum Cell Dimensions 
‘ Coarse 15 levels 58k Ax = Ay=1.7e-3C 
* Refinement studies on grids up to 14.7M points Medium 76 levels B0k Ax = Ay=8.4e-4C 
* Compare with CFL3D, SA model with no circulation correction Fine 17 levels 133 k Ax = Ay =4.2e-4C 
* Mesh convergence sensitive to far-field boundary placement and LE & TE spacings Reference solution 14.7M 1e-7Cx1.25e-5C 
43 “4 
NACA 0012 NACA 0012 
Modified NACA 0012 - Surface pressure mesh convergence Modified NACA 0012 - Skin friction & y* distribution comparison 
pee T'' CEL3D (14.7M cells) 
6 — 15 level (coarse) 
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NACA 0012 


Modified NACA 0012 - Skin friction & y* distribution comparison 


— 15 level (coarse) 
— 16 level (medium) 
— 17 level (fine) 


»CFL3D (14.7M cells) 7 
15 level (coarse) 
16 level (medium) 
17 level (fine) 1500 
Pe : 
500, 


0 02 


04 06 0s 1 0 02 0a 06 08 1 
2/C x/C 

¢ Skin friction shows more sensitivity to resolution 

« Finer meshes do better job near leading edge where physics is most under-resolved 


* HLLC flux responsible for noise — disappears with refinement 
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NACA 0012 


Modified NACA 0012 - Wake surveys of velocity and eddy viscosity 


2/C = 0.999 2/C = 1.001 


+ Experiment and simulation data exist for profiles at 7/C' = 0.999 and 1.001 


+ Examine data at 7/C = 1.001, (similar results at x/C = 0.999) 


NACA 0012 


Wake Velocity Profile 


Wake Eddy Viscosity Profile 


— CFL3D (14.7M cells) 


15 level (coarse) 
16 level (medium) 


— CFL3D (14.7M cells) 
15 level (coarse) 
16 level (medium) 
— 17 level (fine) 


— 17 level (fine) 


NACA 0012 


Modified NACA 0012 - Compare bvp4 with analytic SA-BVP & SA wall function 


Streamwise momentum balance 


— dp/dx 
— udu/dx 
— vdu/dy 
— rhs (sum) 


Skin Friction 
— CFL3D (14.7M cells) (~ 


0.04 


. 2 Ou 
SA-BVP: 2 (w+ uy) 


O05 O05 0 or 02 0 300 0 0 300 000 a Qu) _ Op du dw 
a-velocity * pt ” - F byp4: 5, (+ mo5t) = 92 *? [oe # al 
« Low dissipation inviscid flux helps resolution in cusp - good mesh convergence behavior 4 
i i 0 0.2 04 0.6 08 1 0 0.05 0.1 0.15 
« Lower surface eddy viscosity shows good mesh convergence 2/C «/C 
* Upper surface eddy viscosity not yet mesh converged - literature shows slow convergence (Terms measured at forcing point) 
49 50 
NACA 0012 NACA 0012 


Modified NACA 0012 - Compare bvp4 with analytic SA-BVP & SA wall function 


Streamwise momentum balance 


Skin Friction 


— dpldx 
— uduldx 
— vdu/dy 
— ths (sum) 


0 0.05 01 0.15 
x/IC 


* SA-wall function lags both bvp4 model and SA-BVP skin friction 
¢ SA-BVP does a good job, but misses peak and pressure-side near leading edge 


* Note: RHS relatively small: p: counterbalanced by (uw: + uty) 


Modified NACA 0012 - Mesh refinement near leading edge 


ons 


— cfi3d (14.7M cells) 


EEE 


aie 


= — 17 level (fine) 
— adapted (15 + 2 levels ref.) 


* Mesh refinement important tool for controlling cell counts 
15 level mesh + 2 refinements near leading edge ~ 88 k cells 
* Crshows glitches near refinement interfaces at wall — haven’t dealt with interfaces yet 


NACA 4412 


Modified NACA 4412 with trailing edge separation bubble 


Mo = 0.09 

Rey = 1.52x10° 

Qo = 13.87° 
at 
7 
at 
a 


NACA 4412 


Modified NACA 4412 with trailing edge separation bubble 


Me = 0.09 

Ret = 1.52x10° 

Goo = 13.87° 
Pcl 
-_ 
= 
ee 


* Validation example in the “Extended cases” section of NASA TMR 

* Smooth-body separation bubble near maximum lift conditions 

* Experiment by Coles & Wadcock (1979) with hot-wire velocity profiles 
* Reference data form CFL3D on 897 x 257 grid (230k points) 


NACA 4412 


Modified NACA 4412 with trailing edge separation bubble 


NACA 4412 


Modified NACA 4412 with trailing edge separation bubble 


— 
—_ 
ad 
—_ 
+ 
* Multilevel Cartesian mesh with ~59k cells * Multilevel Cartesian mesh with ~59k cells 
* 1-level of mesh refinement near leading edge * 1-level of mesh refinement near leading edge 
* Leading edge, Ax = 0.1%C, trailing edge Av = 0.2%C, ~1200 cut cells * Leading edge, Ax = 0.1%C, trailing edge Av = 0.2%C, ~1200 cut cells 
56 56 
NACA 4412 NACA 4412 
Modified NACA 4412 — Surface pressure comparison Modified NACA 4412 - Skin friction comparison 
8, oostl —— CFL3D (897 x 257) 
—— Bvp4 
6 — CFL3D (897 x 257) 0.04 —— SA wall function 
5 — Bvp4 
Moo = 0.09 —— SA wall function 
Rez, = 1.52x10° = = ie 
Qo = 13.87° = 0.02 
—_ 
ii 001 
Fadl 
cal 0 
7 0 02 04 0.6 08 
xlC 4 
: _ . * Leading edge very under resolved, y~ ~ 500 
* Good comparison of surface pressure coefficient with both models i 
: Beans nl28, * bvp4 substantially outperforms wall function near leading edge with thin boundary-layer 
* SA-wall function & bvp4 nearly indistinguishable from CFL8D results & steep pressure gradient 
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Modified NACA 4412 - Skin friction comparison Modified NACA 4412 - Velocity comparison near separation 
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* bvp4 substantially outperforms wall function near leading edge with thin boundary- 
& steep pressure gradient 
* bvp4 predicts separation location within 1% of mesh resolved CFL3D resolved 


* Noise in bvp4 due to interpolation of (wv ux & v uy) at forcing point 


ayer 


Locations of hot-wire surveys in experiment (Coles & Wadcock, 1979) 
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Modified NACA 4412 - Velocity comparison near separation 
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~~~ CFL3D (dashed) == CFLAD (dashed): 


0.09] — BVP4 (solid) 


0.08 


x= 0.9528 
x = 0.8973 
x = 0.8418 
x = 0.7863 


0.07 
0.06 


04 x= 0.6573 


Vertical distance from wall 
Vertical distance from wall 


-~002 0 0.02 0.04 0.06 0.08 0.1 -0.012 -0.008 -0.004 a) 0.004 
X-—velocity y-—velocity 


* Good prediction of both x and y components of velocity through separation bubble 
* Vertical velocity about an order of magnitude smaller than horizontal 


* Slight “viscous overshoot” due to coarseness of Cartesian mesh, Ax = Ay = 0.2%C 
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Summary 


Atmospheric propagation and ground effects modeling 
+ presented V&V studies for a new ODE-based wall model for RANS equations 


+ demonstrated for several well-studied flows including smooth body separation 


* bvp4 model: 


+ Solves coupled solves a coupled set of ODEs posed as two-point boundary value problems for 
the streamwise velocity and the turbulent viscosity 


+ includes both the streamwise pressure gradient and the momentum balance valid farther from the 
wall 


permits wall spacing on the Cartesian mesh 4 to 8 x coarser than with analytic wall functions, 
order of magnitude farther out than analytic wall functions 
+ wall model itself about 2-3x the computational cost of analytic wf’s on same mesh 


+ Can be applied in body-fitted or non-body fitted meshes ¥ 


TODO 


+ Paper presents counting arguments for 3D 


Questions’ 


NACA 4412 


Modified NACA 4412 - SA model vs experiment data near separation 
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* Very sensitive example - near CL max 

« SA makes some error in velocity profiles due to mis-prediction of separation location 
* Different turbulence models show up to 10%C variation in separation location 

« Vertical velocity even more sensitive, since magnitudes ~10 x smaller 


NACA 0012 


Modified NACA 0012 - Integrated force coefficients 


+ TMR has loads data available for several codes 


* Excellent prediction of viscous drag. Even 
coarsest mesh (15lev) is within 2 counts 


* Net lift and drag not as good due to inviscid 
regions of flow 


* TMR documents sensitivity to far field boundary 
(100C' vs 500C for reference) 


aif . CaViscous Ca Pressure CaTotal C1 Total 
* Trailing edge spacing 


too large: ee 0.00621 | 0.00607 | 0.01227 | 1.0908 
.7M point reference 
~500 x coarser than 


17 Lev 
reference on adapted uniform wall pacen, 1aakeets| 0-20611 | 0.00767 | 0.01378 | 1.1120 
grid Adapted 


0.00607 0.00751 0.01358 1.1416 


15 lev + 2 near LE, 88k cells. 


t Data from CFL3D with SA model on “family II” grid, no point vortex correction & 2nd-order 
turbulent advection. 


